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This  report  summarizes  the  work  conducted  during  the  period  of  Jan.  1,  2002  to  March 
31,  2005.  During  the  grant-supported  period,  the  PI  continues  to  work  on  modeling  and 
computation  of  flows  of  liquid  crystal  polymers  (LCPs)  and  nanocomposites  guided  by 
Air  Force  contact  Dr.  Richard  Vaia  of  material  directorate  WPAFB.  The  activities 
include  developing  kinetic  theories  for  uniaxial  molecules  with  finite  extensibility,  flows 
of  cholesterics,  nonhomogeneous  flows  accounting  for  stress  induced  polymer  migration, 
as  well  as  theories  for  biaxial  liquid  crystals  and  nano-particle  suspension  flows.  A 
continuum  theory  for  flow  of  cholesteric  liquid  crystals  is  also  derived  to  explore  the 
mesostructure  in  weak  flows.  Works  on  nematic  rigid  polymers  are  refined  with  the 
kinetic  theory  calculations.  Materials  properties  for  nanocomposites  are  explored  through 
homogenization  techniques  using  our  mesoscopic  structure  calculations  and  an  effective 
medium  theory.  The  main  collaborators  in  this  period  are  Prof.  M.  G.  Forest  of  UNC 
Chapel  Hill,  Prof.  Ruhai  Zhou  (former  Postdoc  at  UNC  and  now  at  Old  Dominion),  Prof. 
Bechtel  of  Ohio  State  University,  Professor  M  C  Calderer  of  University  of  Minnesota, 
Prof.  Hong  Zhou  of  Naval  Post  Graduate  School,  Prof.  Rob  lipton  of  LSU,  and  my 
graduate  students  Zhenlu  Cui  and  Sarthok  Sircar  of  FSU.  (Zhenlu  Cui  has  graduated  this 
summer  with  a  Ph.  D.  degree  in  mathematics.)  In  the  following,  I  summarize  the 
completed  work  and  the  work-in-progress  briefly  and  list  the  publications  within  the 
period.  For  details,  please  refer  to  our  publications. 

Multiscale  modeling  and  study  of  flow-orientation  coupling  in  the  shear 
flow  of  liquid  crystal  polymers. 

The  Pi  and  his  collaborators  (M.  G.  Forest  and  R.  Zhou)  continue  refined  analysis  on  the 
extended  Doi-Hess  kinetic  theory  for  mondomain  shear  flows  and  general  planar  linear 
flows  of  liquid  crystal  polymers.  The  aim  here  is  to  provide  a  comprehensive 
understanding  of  the  monodomain  flow  behavior  governed  by  the  kinetic  theory.  In  a 
series  of  papers,  we  thoroughly  examined  the  monodomain  dynamics  and  characterized 
completely  the  complex  dynamics  governed  by  the  kinetic  theory  using  a  Galerkin 
method  based  on  the  spherical  harmonic  expansions  and  a  mesoscopic  closure  model.  In 
the  dilute  regime,  we  solved  the  governing  equation  asymptotically  in  the  weak  shear 
limit  and  give  the  explicit  scaling  behavior  of  the  solution  with  respect  to  the  important 
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flow  parameters  as  well  as  the  geometric  parameters  of  the  molecules.  The 
comprehensive  phase  diagram  obtained  using  both  closure  and  kinetic  theory  provides  an 
atlas  for  further  studies  of  the  nonhomogeneous  flow  behavior  of  LCPs.  Following  the 
monodomain  road  map,  we  have  identified  a  variety  of  temporal-spatial  structures  and 
patterns  in  sheared  nematic  LCPs  of  nonhomogeneous  flows.  The  newly  identified 
chaotic  behavior  in  monodomain  is  also  observed  in  the  nonhomogeneous  flow  through 
both  temporal  and  spatial  ramifications. 

In  additional  to  the  sheared  nematic  LCPs,  we  also  examined  the  kinetic  theory  in  general 
linear  flows  and  discovered  that  the  flow  map  can  be  completely  recovered  from  our 
results  for  sheared  nematics.  The  crucial  link  is  the  extended  geometry  parameter  a.  In 
other  words,  the  general  linear  flow,  including  the  important  4-roll-mills,  can  be  obtained 
by  rescaling  the  results  for  the  sheared  nematics.  In  two  papers,  we  explored  the  basins 
of  attractors  and  the  chaotic  boundaries  using  kinetic  and  closure  models,  respectively. 
The  results  are  published  in  PRL  and  Rheological  Acta,  respectively. 

Our  current  activities  in  this  direction  are  (i).  summarize  our  extensive  results  with 
respect  to  the  nonhomogeneous  mesoscopic  structure  and  have  them  published  and  to 
compute  the  governing  equation  obtained  via  closures  in  multidimensional  geometry 
aiming  at  understanding  the  flow-orientation  interaction  in  sheared  nonhomogeneous 
nematics  and  how  the  temporal  structure  in  the  monodomain  interact  with  the  spatial 
degree  of  freedom.  Ref  [3-6, 9-1 1,13,18,20,24,25,33]. 

Multiscale  models  for  uniaxial  molecules  of  finite  extensibility  and 
nonhomogeneous  liquid  crystal  polymers. 

The  kinetic  theory  developed  in  (Wang,  2002)  has  been  extended  to  account  for  the 
spatial  inhomogeneity  of  the  lcp  by  allowing  the  density  of  the  material  to  vary  in  space 
(Q  Wang,  M  G  Forest,  and  R  Zhou,  2004).  This  additional  degree  of  freedom  is  amount 
to  introduce  another  order  parameter  for  the  spatial  distribution  of  the  center  of  mass  of 
the  liquid  crystal  polymers.  Some  spatial  patterns  related  to  the  center  of  mass  variation 
are  therefore  plausible.  Indeed,  we  can  solve  the  governing  equations  in  the  moment- 
averaged  approximation  to  show  that  a  periodic  density  pattern  may  form  while  the 
averaged  molecular  orientation  is  attained  absent  of  flows.  The  potential  layering  is 
naturally  built  into  the  kinetic  theory.  We  remark  that  the  density  variation  of  LCPs  is  not 
modeled  in  the  classical  theory  of  Ericksen-Leslie  or  any  other  continuum  mechanical 
theories  nor  in  the  phenomenological  theories  of  Tsuji-Rey's  etc.  In  our  most  recent 
studies,  we  found  the  spatial  variation  of  densities  and  the  associated  translational 
diffusion  is  intrinsically  related  to  the  long-range  elastic  distortion  of  the  molecular 
orientation. 

Lately,  the  theory  is  further  refined  to  include  the  stress  induced  polymer  migration,  a 
phenomenon  very  important  to  the  translational  diffusion  of  polymer  molecules  in 
nonhomogeneous  flows,  especially,  near  the  wall.  Detailed  analyses  on  the  mathematical 
structure  of  the  model  as  well  as  the  rheological  consequence  of  the  effect  are  currently 
summarized  in  a  preprint. 
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An  approximate  theory  for  uniaxial  molecules  is  being  examined  with  appropriate 
parameters  to  evaluate  the  effect  of  molecular  stretching  on  the  rheological  behavior  of 
the  flowing  material.  Experimental  rheologists  has  long  suspected  that  the  conformational 
change  of  LCP  molecules  in  the  flow  may  have  been  one  of  the  key  ingredients  that  has 
been  missing  in  the  rigid  rod  or  disk  models  for  LCPs.  With  this  new  model,  we  can 
capture  the  three  distinctive  viscosity  regime,  eluded  the  previous  lcp  theories.  We  hope 
we  will  be  able  to  give  more  detailed  explanations  for  this  speculation  soon.  Ref 
[1,2,8,12,22,28], 

Material  property  characterization  through  mesoscopic  structure 
calculations 

Through  interaction  with  Dr  Vaia  at  WPAFB  and  joined  by  Prof.  Lipton  of  LSU,  we 
carried  out  a  series  study  on  the  conductivity  property  of  polymer-particle-inclusion 
nanocomposites.  We  employed  a  homogenization  technique  to  arrive  at  the  effective 
conductivity  property  of  the  nanocomposites  described  by  the  mesoscopic  structure 
tensor  of  the  nanocomposites.  The  mesoscopic  structure  is  the  output  of  our 
nonhomogeneous  theory  for  nanocomposites.  The  conductivity  properties  at  various 
material  parameter  regimes  are  investigated.  Spatially  inhomogeneous  morphology  of  the 
nanocomposites  is  explored  recently  giving  rise  to  an  effective  medium  property 
throughout  the  physical  domain  where  the  medium  resides.  Ref  [  1 9,2 1 ,23 ,27] . 

Modeling  incompressible  fluids  through  continuum  mechanics  principles 

This  is  a  joint  work  with  Prof.  Stephen  Bechtel  of  Ohio  State  and  Frank  Rooney.  In  a 
continuing  effort  to  model  incompressible  fluids  with  prescribed  density  variations,  we 
examined  the  models  we  developed  for  incompressible  fluids  with  various  prescribed 
density  variations  through  mechanics  principles.  We  clarified  the  definition  of  the 
pressure  for  these  types  of  constraint  theories  and  studied  the  flow  and  temperature 
profiles  in  plane  Couette  and  Poiseuille  flows.  The  results  are  published  in  International 
Journal  of  Engineering  Science  and  Physics  of  Fluids.  Ref  [14,  17]. 

Multiscale  kinetic  theories  for  flows  of  Biaxial  liquid  crystal  polymers 

Given  the  rising  interests  in  the  modeling  of  nanofluids  of  biaxial  constituents  such  as  the 
flow  of  liquid  crystal  polymers,  we  developed  a  hydrodynamic  theory  for  the  biaxial 
liquid  crystal  polymer  by  modeling  it  as  an  ellipsoid  suspended  in  viscous  solvent.  The 
excluded  volume  interaction  potential  is  calculated  from  the  first  principle.  Detailed 
analysis  is  underway  and  will  be  reported  in  the  next  annual  report  to  AFOSR. 


Permeation  mode  in  weak  flows  of  cholesteric  liquid  crystal  polymers 
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We  have  obtained  the  mesoscopic  structure  using  a  coarse-grained  approximation  and 
homogenization  technique  for  very  small  pitch  cholesterics.  We  obtained  the  order 
parameter  variation  across  the  flow  field  as  well  as  their  influence  to  the  major  director 
orientation  in  the  flow.  The  work  has  been  submitted  to  DCDS-series  B  and  J.  of  Non- 
Newtonian  Fluid  Mechanics  (2005).  Ref  [28,29]. 

On-going  projects 

We  are  in  the  process  of  preparing  submissions  on  our  work  about  the  monodomain 
dynamical  study  of  models  for  lcp  molecules  of  finite  extensibility  and  mesostructure 
simulation  of  nonhomogeneous  lcps  with  density  variations.  Our  3-D  code  for  the 
simulation  of  shear  nematic  rigid  polymers  is  in  the  process  to  be  parallelized  and  run  on 
super-computers. 

The  research  activities  on  the  projects  funded  by  the  award  are  proceeding  smoothly 
according  to  the  plan.  Additional  projects  closely  related  to  the  original  proposed  are 
identified  and  added  to  our  research  activities.  We  anticipate  the  slightly  expanded  scope 
of  research  will  facilitate  our  efforts  for  the  overall  objective  as  well  as  meet  the  need  of 
the  Air  Force.  We  will  report  out  progress  on  the  on-going  projects  in  our  next  technical 
report. 


Publications  during  the  period: 


1 .  Q.  Wang,  "A  hydrodynamic  theory  of  nematic  liquid  crystalline  polymers 
of  different  configurations",  Journal  of  Chemical  Physics,  116  (2002)  pp. 
9120-9136, 

2.  Qi  Wang,  Weinan  E,  Chun  Liu,  and  Pingwen  Zhang,  "  Kinetic  theories  for 
flows  of  nonhomogeneous  rodlike  liquid  crystalline  polymers  with  a 
nonlocal  intermolecular  potential",  Physical  Review  E,  65(5)  (2002),  pp. 
0515041-0515047. 

3.  M.  G.  Forest  and  Q.  Wang,  "  Monodomain  response  of  finite-aspect-ratio 
macromolecules  in  shear  and  related  linear  flows",  Rheological  Acta  42 
(2003),  20-46. 

4.  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  ''Symmetries  of  the  Doi  kinetic 
theory  for  nematic  polymers  of  finite  and  infinite  aspect  ratio:  at  rest  and 
in  linear  flows,"  Physical  Review  E  66(3)  (2003),  P031712. 

5.  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  "Explicit  flow-aligned  orientational 
distribution  function  for  dilute  nematic  polymers  in  weak  shear," 
RHEOLOGY  AND  FLUID  MECHANICS  OF  NONLINEAR  MATERIALS, 
edited  by  D.  A.  Siginer,  D.  DeKee,  and  S.  Bakhatiyarov,  ASME,  New 
York,  2002. 

6.  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  "Full  tensor  alignment  criteria  for 
sheared  nematic  polymers",  J.  of  Rheology  47(1)  (2003),  pp.  105-127. 
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7.  S.  E.  Bechtel,  M.  G.  Forest,  F.  J.  Rooney,  Q.  Wang,  "Thermal  Expansion 
Models  of  Viscous  Fluids  Based  on  Limits  of  Free  Energy",  Physics  of 
Fluids,  15(9)  (2003),  00.  2681-2693. 

8.  Q.  Wang,  M.  G.  Forest  and  R.  Zhou,  ”A  hydrodynamic  theory  for  solutions 

of  nonhomogeneous  nematic  liquid  crystalline  polymers  with  density 
variations,"  J.  of  Fluid  Engineering,  126  (2004),  pp180-188. 

9.  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  "The  weak  shear  phase  diagram  for 
nematic  polymers,"  Rheological  Acta,  43(1)  (2004),  pp.  17-37. 

10.  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  "Scaling  behavior  of  kinetic 
orientational  distributions  for  dilute  nematic  polymers  in  weak  shear,"  J. 
Non-Newtonian  Fluid  Mechanics,  116  (2004),  pp.  183-204. 

11.  M.  G.  Forest,  Q.  Wang,  H.  Zhou,  R.  Zhou,  "Scaling  Structure  scaling 
properties  of  confined  nematic  polymers  in  plane  Couette  cells:  the  weak 
flow  limit,"  J.  Rheology,  48(1)  (2004),  pp.  175-1 92. 

12.  M.  C.  Calderer,  M.  G.  Forest,  and  Q.  Wang,  "Kinetic  Theories  and 
Mesoscopic  Models  for  Solutions  of  Nonhomogeneous  Liquid  Crystal 
Polymers,"  J.  Non-Newtonain  Fluid  Mechanics,  120  (2004),  pp.  69-78. 

13.  M.  G.  Forest,  Q.  Wang,  R.  Zhou,  and  E.  Coate,  "Monodomain  response  of 
arbitrary  aspect  ratio  nematic  polymers  in  general  linear  planar  flows,"  J 
of  Non-Newtonian  Fluid  Mechanics,  118(1)  (2004),  pp.  17-31. 

14.  S.  E.  Bechtel,  F.  Rooney,  and  Q.  Wang,  "A  thermodynamic  definition  of 
pressure  for  incompressible  viscous  fluids,"  International  J.  of 
Engineering  Science,  (42)  19-20  (2004),  pp.  1987-1994. 

15.  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  The  flow-phase  diagram  of  Doi- 
Hess  theory  for  sheared  nematic  polymers  II:  finite  shear  rates, 
Rheological  Acta,  Vol  44(1)  (2004),  pp.  80-93. 

16.  X.  Zheng,  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  "Likelohood  and  expected 
-time  statistica  of  monodomain  attractors  in  sheared  discotic  and  rodlike 
nematic  polymers,"  Rheololical  Acta,  Vol  44(3)  (2005),  pp.  219-234. 

17.  S.  E.  Bechtel,  M.  G.  Forest,  F.  Rooney,  and  Q.  Wang,  "Investigation  of 
simplified  thermal  expansion  models  compressible  Newtonian  Fluids 
applied  to  nonisothermal  plane  Couette  and  Poisuille  flows,"  Physics  of 
Fluids,  in  press  2004. 

18.  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  "Chaotic  boundaries  of  nematic 
polymers  in  mixed  shear  and  extensional  flows,"  Physical  Review  Letters, 
93  (8)  (2004),  pp.  088301-088305. 

19.  X.  Zheng,  M.  G.  Forest,  R.  Lipton,  R  Zhou,  Q.  Wang,  “Exact  scaling  laws 
for  electrical  conductivity  properties  of  nematic  polymer  nano-composite 
monodomains,”  Advanced  Functional  Materials,  in  press,  2004. 

20.  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  “Kinetic  structure  simulations  of 
nematic  polymers  in  plane  Couette  cells,  I:  The  algorithm  and 
benchmarks,”  SIAM  MMS,  in  press,  2004. 

21.  H.  Zhou,  M.  G.  Forest,  X.  Zheng,  Q.  Wang,  R.  Lipton,  “Extension- 
enhanced  conductivity  of  liquid  crystalline  polymer  nano-composites,” 
Macromolecular  Symposia,  in  press,  2004. 
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22.  M.  G.  Forest  and  Q.  Wang,  Hydrodynamic  theories  for  blends  of  flexible 

polymer  and  nematic  polymers,  Physical  Review  E,  in  press  2004. 

23.  M.  G.  Forest,  R.  Zhou,  Qi  Wang ,  X.  Zheng,  R.  Lipton,  “Anisotropy  and 
Heterogeneity  of  Nematic  Polymer  Nano-Composite  Film  Properties,”  to 
appear  in  IMA  volume  on  Soft  Matter,  2005. 

24.  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  Kinetic  structure  simulations  of 

nematic  polymers  in  plane  Couette  cells,  II.  to  appear  SIAM  MMS,  2005. 

25.  Z.  Cui,  M.  G.  Forest,  and  Q.  Wang,  "Effects  of  the  long-range  anisotropic 
elasticity  in  weak  plane  Couette  and  Poisuille  flows  of  nematic  liquid 
crystal  polymers,"  to  appear  Siam  Journal  on  Applied  Math,  2005. 

26.  P.  B.  Chase,  A.  Kataoka,  B.  C.  W.  Tanner,  J.  M.  Macpherson,  X.  Xu,  Q. 
Wang,  M.  Reginier  and  T.  Daniel,  “Spatially  explicit,  nanomechanical 
models  of  the  muscle  half  sarcomere:  Implications  for  mechanical  tuning 
in  atrophy  and  fatigue,”  to  appear  in  Astraunica  Acta,  2005. 

27.  M.  G.  Forest,  R.  Zhou,  Qi  Wang ,  X.  Zheng,  R.  Lipton,  “Anisotropy  and 

dynamics  ranges  in  effective  properties  of  sheared  nematic  polymer 
nano-composites,”  Advanced  Functional  Materials,  in  press  2005. 

PAPERS  SUBMITTED 


28.  Z.  Cui  and  Q.  Wang,  "A  continuum  mechanics  model  for  flows  of  chiral 
nematic  polymers  and  permeation  flows,"  submitted  to  J.  of  Non- 
Newtonian  Fluid  Mechanics,  2005. 

29.  Z.  Cui,  M.  C.  Calderer,  Q.  Wang,  "A  kinetic  theory  for  flows  of  cholesteric 
liquid  crystal  polymers,"  submitted  to  Discrete  and  Continuous  Dynamical 
Systems-Series  B,  2005. 

30.  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  “Extending  the  magnetic/electric 
field-hydrodynamic  analogy  of  rod  &  platelet  suspensions  to  strongly 
coupled,”  submitted  to  Phys.  Rev.  Lett.  2005. 

31.  H.  Zhou,  H.  Wang,  M.  G.  forest,  and  Q.  Wang,  “A  new  proof  on  uniaxial 
equilibria  of  Smoluchowski  equation  for  rodlike  nematic  polymers,” 
submitted  to  Nonlinearity,  2005. 

32.  Q.  Wang,  S.  Saircar,  and  H.  Zhou,  “Solutions  of  the  Smoluchowski 
equation  for  nematic  polymers  under  imposed  fields,”  submitted  to 
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